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Abstract 

The hystcrctic behavior of many-particle systems with non-convex free energy can be modeled 
by nonlocal Fokker-Planck equations that involve two small parameters and are driven by a time- 
dependent constraint. In this paper we consider the fast reaction regime related to Kramers- type 
phase transitions and prove that the dynamics in the small-parameter limit can be described by a 
^ rate-independent evolution equation. To this end we derive mass-dissipation estimates from Muck- 

^ enhoupt constants, establish dynamical peak-stability estimates, and employ moment estimates 

Q that encode large deviations results. 
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1 Introduction 



It is an ubiquitous and intriguing question in the mathematical analysis under which conditions 
the dynamics of a given high-dimensional systems with small parameters can be described by low- 
dimcnsional, reduced evolution equations. In this paper we answer this question, at least partially, 
for a particular example, namely the Fokker-Planck equation 

TdtQ{t, x) = d^(u^d^Q{t, x) + {H'{x) - ait)) git, x)), (FPi) 

where r and u are the small parameters and x G M is a one-dimensional state variable. Moreover, 
H is supposed to be a double-well potential and cr is a dynamical multiplier chosen such that the 
solution complies with 

[ XQit,x)dx = iit), (FP2) 
Jr 

where £ is a prescribed control function. This dynamical constraint is, for admissible initial data, 
equivalent to the mean-field formula 

ait)= [ H'ix)Qit,x)dx + Tlit), (FP^) 

which turns (FPi) into a nonlocal, nonlinear, and non-autonomous PDE. 

Nonlocal Fokkcr-Planck equations like (FPi) + (FP2) have been introduced in [DGHll] in order to 
model the hysteretic behavior of many-particle storage systems such as modern Lithium-ion batteries 
(for the physical background, we also refer to [DJG+10]). In this context, x G M describes the 
thermodynamic state of a single particle (nano-particle made of iron- phosphate in the battery case) , 
H is the free energy of each particle, and i' accounts for entropic effects. Moreover, g is the probability 
density of a many-particle ensemble and the dynamical control i reflects that the whole system is 
driven by some external process (charging or discharging of the battery). 

Since H is non-convex, the dynamics of (FPi)-|-(FP2) can be rather involved as they are related 
to three different time scales, namely the small relaxation time r, the time scale of the control i 
(which is supposed to be of order 1), and the Kramers scale rexp (/i((t)/z/^), which corresponds 
to probabilistic transitions between the different wells of a time-dependent effective potential with 
energy barrier hia). The different dynamical regimes for < i/,t have been investigated by the 
authors in [HNV12] using formal asymptotic analysis. 

In this paper we restrict our considerations to the fast reaction regime, that means we suppose 
< ^ 1 and assume that r is coupled to by a certain exponential scaling law implying < t <^ u. 
In the most simple and prototypical case, this scaling law reads 




where is some given parameter that is positive but smaller than a certain threshold /ithres- We 
emphasize that there exists also also a slow reaction regime corresponding toO<i^^r^l, 
but then the dynamics is more complicated and neither related to rate-independent evolution nor 
Kramers-type phase transitions, see the discussion in [HNV12]. 

Our main result is the proof that the microscopic PDE (FPi)-|-(FP2) can be replaced, as ^ 0, by 
a low-dimensional dynamical system, which turns out to be rate-independent and exhibits hysteresis. 
These macroscopic equations govern the evolution of the multiplier a and the phase fraction jj,, which 
is defined by 

fiit) = / git, x)dx — / git, x) dx, 

J right stable region J left stable region 

where 'stable region' refers to a connected component of {x : H"ix) > 0}. 
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The micro-to-macro transition studied here is similar to those in [PT05, A/[iellb, MT12], which 
Hkewise derive macroscopic models for hysteric behaviour from microscopic gradient flows with non- 
convex energy and external driving. Our microscopic system, however, is different as it involves 
the diffusive term v^d^g, which causes specific effects and necessitates the use of different methods. 
More precisely, the dominant effect in the fast reaction regime of nonlocal Fokker-Planck equations 
are Kramers-type phase transitions, which describe that particles can pass through the spinodal 
region {x : H"{x) < 0} due to stochastic fluctuations. 

The key observation in our context is that Kramers-type phase transitions can manifest on the 
macroscopic scale only if the dynamical multiplier a attains one of two critical values and o"*, 
which are completely determined by H and h^, because otherwise the corresponding microscopic 
mass flux is either too small or too large. The limit dynamics for — > is therefore completely 
characterized by the flow rule 

fi{t) < for a{t) = a^, fi{t) > for a{t) = cr*, fi{t) = otherwise, 

and pointwise relations C{i{t), (j{t), ^(i)) = that encode the dynamical constraint. These findings 
can be summarized as follows. 

Main result. Under natural assumptions on H , the control I, and the initial data, the triple {i, a, ji) 
satisfies in the limit — >■ a closed rate-independent evolution equation with hysteresis. Moreover, 
the limit solution is unique provided that the initial data are well-prepared. 

The rest of the paper is organized as follows. In §2 we give a more detailed introduction into 
the problem. In particular, in §2.1 and §2.2 we specify our assumptions and review the existence 
theory for (FPi)-|-(FP2) with arbitrary z/, r > as it is developed in Appendix A. Moreover, in §2.3 
we heuristically explain the key dynamical features in the fast reaction regime and proceed with a 
precise formulation of the limit model in §2.4. 

A major part of our analytical work is contained in §3. Specifically, we establish mass-dissipation 
estimates in §3.1 and derive in §3.2 conditional results for the dynamical stability of localized peaks. 
Afterwards we study the mass transfer between the two stable regions in §3.3 and §3.4. 

In §4 we pass to the limit z/ — > 0. We continue our investigations concerning the dynamical 
stability of peaks in §4.1 and obtain uniform Lipschitz estimates for the multiplier a in §4.2. These 
ingredients finally enable us to prove our main result in §4.3, see Theorems 29 and 30. 

2 Preliminaries 

In this section we introduce our assumption on H, i, and the initial data, and summarize some 
important properties of solutions to the non-local Fokker-Plank equation. Moreover, we discuss the 
dynamics in the fast reaction regime on a heuristic level and formulate the rate-independent limit 
model. 

2.1 Assumptions on the potential 

Throughout this paper we assume that H \s & double- well potential with the following properties, 
see Figure 1 for an illustration. 

Assumption 1 (properties of H). 

1. H is three times continuously differentiable, attains a local maximum at x = and the global 
minimum at precisely two points. 

2. H" has only two zeros x*, x* with x* < < x* ; we set a* = H'{x*) and a* = H'{x^) and this 
implies a* < < a* . 

3. H' is asymptotically linear in the sense of \imx-^±^ H'" (x) = and limx-:f±oo H" (x) = c± for 
some constants c± . 
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Figure 1: Example of a double- well potential H that satisfies assumption 1 with cr# and cr'^ as in Assumption 
4. The shaded regions illustrate the spinodal (or unstable) interval {x*, a;*). 



The assumption that the two wells of H are global minima is not crucial and can always be 
guaranteed by means of elementary transformations. In fact, (FPi) and (FP2) are, for any given 
c G M, invariant under H H + cx, a a + c. Moreover, by an appropriate shift in x we can always 
ensure that the local maximum is attained at x = 0. The assumption that H grows quadratically at 
infinity is of course more restrictive and made in order to keep the presentation as simple as possible. 
We expect, however, that our convergence result is also true for more general double- well potentials 
H provided that these grow superquadratically or that the initial data decay sufficiently fast. 



As a direct consequence of Assumption 1 we can introduce three functions X_, Xq, and X+ such 
that H' o Xj = id. 

Remark 2 (functions X_, Xq, and X^). The inverse of H' has three strictly monotone and differ- 
entiable branches 

X_ : [—00, a*) — >■ (—00, X*] , Xq : [cr*, a*] [x*, x*] , X+ : [cr*, +00) — >■ [x*, +00) . 
In particular, we have 

1. X+{a) — X_((t) > c for a all a & [a*, a*], 

2. c< X'^{a) < Ce for all a G [ct* + e, a* - e], 

3. \g2 — (y\\< C\X±{(j-2) — X-|-((7i)| for all o\,G2 in the domain of X±, 
for any e with < e < \{(t* — ct*) and some constants c, C and Ce. 

In order to describe Kramers-type phase transitions, we further introduce the effective potential 

H^{x) := H{x) — ax, 

and define two functions /i+ : (cr*, cr*) ^ M by 

h±{a) := H^iXoia)) - H„{X±{a)). 

These definitions are motivated by the many-particle interpretation of (FPi). In fact, for frozen a 
the particles diffuse in the effective potential and the energy barriers h- and appear explicitly in 
Kramer's formula for the mass fluxes between the two wells of H^, see Figure 2 and the discussion 
in §2.3. 

Remark 3 (properties oih±). The functions h- and h+ are well-defined and smooth on the interval 
[a^, a*] with /i-(0) = /i+(0) > 0. Moreover, h- is strictly decreasing with h-{a*) = and hj^ is 
strictly increasing with h-^{a^) = 0. 

We finally describe the coupling between r and v and introduce the values (T# and cr#. 
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Figure 2: Cartoons of the effective potential H„ with a* < a < and the functions h- (sohd hne) and /i+ 
(dashed line). The values cr# and tr* are defined by /i_ (c*) = /i+(cr#) = with h^ = — liniy^o ^"^ ^'o.t. 

Assumption 4 (coupling between r and v). The parameter r is positive, depends on v, and satisfies 

2 1 ''^O , 

v mr > — /i# 

for some with < h^ < ^thres •= ^±(0)- In particular, there exist and a"^ such that 

<a# <0 <a* < a*, h# = h-{a*) = h+(a^), 
and hence < /ithres < /i+(<7'^)}. 
2.2 Existence and properties of solutions 

It is well established, see [JK097, JK098], that the linear Fokker-Planck equation without dynamical 
constraint - that is (FPi) with a{t) = - is the Wasserstein gradient flow to the energy 



£(t) := / ^(t, a;)ln^(i, dx + / H{x)Q{t, x)dx. 
J Jr 



(1) 



Similarly, the non-driven variant of the nonlocal Fokker-Planck equations - that is (FPi)+(FP2) 
with £{t) = - can be regarded as the Wasserstein gradient flows for £ restricted to the constraint 
manifold J^gdx = £, and we easily verify that the corresponding dissipation is given by 



m ■■= I 



(u^d^git, x) + {H'{x) - a{t))Q{t, x))" 
g{t, x) 



dx. 



(2) 



In the general case £ ^ 0, however, the energy is no longer strictly decreasing but satisfies 

T£{t) = -V{t) + Ta{t)i{t). (3) 

In particular, wc have d£ < adi along each trajectory, and this reflects the second law of thermo- 
dynamics for the free energy of the many-particle ensemble in the presence of the dynamical control. 
The energy-dissipation estimate (3) is essential for passing to the limit ^ as it reveals that the 
dissipation V is very small with respect to the L^-norm and hence, loosely speaking, also small at 
most of the times. For linear Fokker-Planck equations without constraint, the underlying gradient 
structure can be used to establish F-convergence as r — )■ 0. The resulting evolution equation is a 
one-dimensional reaction ODE for the phase fraction /x and equivalent to Kramers' celebrated for- 
mula, sec [PSVIO, AMP+11, HNll]. However, it is not clear to us whether this variational approach 
can be adapted to the present case with dynamical constraint; the methods developed here employ 
the estimate for V but make no further use of the gradient flow interpretation of (FPi)-|-(FP2). 



Since the system (FPi)-|-(FP2) is a nonlinear and nonlocal PDE, it is not clear a priori that the 
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initial value problem is well-posed in an appropriate function space. In the case of a bounded 
spatial domain and Neumann boundary conditions, the existence and uniqueness of solutions has 
been established in [Hutl2, DHM+11] using an L^-setting for q with q > 1. Since here we are 
interested in solutions that are defined on the whole real axis, we sketch an alternative existence and 
uniqueness proof in Appendix A. The key idea there is to obtain solutions as unique fixed points of a 
rather natural iteration scheme on the state space of all probability measures with bounded variance. 
Moreover, adapting standard techniques for parabolic PDE we derive several bounds to reveal how 
these solutions depend on u. 

Our assumptions and key findings concerning the existence and regularity of solutions to the 
nonlocal Fokker-Planck equation be can be summarized as follows. 

Assumption 5 (dynamical control I). The final time T with < T < oo is independent of v. The 
control I is also independent of v and twice continuously differentiahle on [0, T]. In particular, we 
have 

sup (\t{t)\ + \t{t)\ + \i{t){) <c 

tG[0,T] ^ ^ 

for some constant C independent of v. 

Assumption 6 (initial data). The initial data are nonnegative and satisfy 

/ giO,x)dx = l, / xq{Q, x)dx = i{0), / x"^ giO, x) dx < C 
Jr Jr Jr 

for some constant C independent of v. 

Lemma 7 (existence and properties of solution). For any v with < < 1 and given initial data 
there exists a unique solution g to the initial value problem (FPi)-f(FP2) which is nonnegative and 
smooth for t > 0, and satisfies 

/ g{t,x)dx = l, / xg{t, x) dx = i{t) 
Jr Jr 

for all t G [0, T] . Moreover, each solution satisfies 

sup (\a{t)\+ [ x'^g{t, x)dx^ + sup \\g{t, ■)\\oo + t'^ [ V{t)<C 

<G[o,T] \ Jr J te[t*,T] Jtt 

with for some constant C which depends only on H, I and J^x'^g{0, x)dx. 

Proof. All claims follow from Proposition 31 and Proposition 32 in Appendix A. □ 

The assertions of Lemma 7 reflect the existence of two small transient time scales. At first we 
have to wait for the time before we can guarantee that \\g(t, ■)\\oo < C/iy"^ and T>{t) dt < Ct. 
The first estimate is needed within §3 in order to show that no mass can penetrate the spinodal region 
from outside, and that there is no mass flux through the spinodal region for subcritical a G (a^, cr*) . 
Furthermore, it is in general not before a time of order t^~^ that the dissipation V{t) is eventually 
smaller than (the exponent < /3 < 1 will be identified below). In §4 we prove that the solutions 
to the nonlocal Fokker-Planck equations behave nicely after the second time, even though we are not 
able to exclude that T>{t) becomes large (again) at some later time. 

The initial transient regime corresponds to very fast relaxation processes during which the system 
dissipates a large amount of energy leading to rapid changes of especially the multiplier a and the 
phase fraction /x. For generic initial data, we therefore expect to find several limit solutions as 
1/ — >■ depending on the microscopic details of the initial data. The only possibility to avoid such 
non-uniqueness is to start with well-prepared initial data. 



6 



Definition 8 (well-prepared initial data). The initial data from Assumption 6 are well-prepared, if 

they additionally satisfy 

A\q{^, •)lloo + r-^2^(0)<C, 
for some constant C independent of v, and if we have 

Cr(0) > (Jini 

for some crini G M. 

Remark 9. For well prepared initial data we can choose t* = in Lemma 7. Moreover, we have 

/n \ "^"^0 ^ 1 ~ A^ini j / \ i -'- Mini j / s 

^(0, X) > Qini := ^ (5x_((Tini)(^) + 2 

in the sense o^ weak* convergence of measures, where Sx denotes the Dirac distribution at X G U. 
and /Ltini := f^^ Qim{x) dx - f^^°° gini{x) dx. 

Proof The assertions follow from Remark 33 and the mass dissipation estimates formulated in 
Lemma 17 and Lemma 18. □ 

2.3 Heuristic description of the fast reaction regime 

In order to highlight the key ideas for our convergence proof, we now give an informal overview on 
the effective dynamics for <C 1. For numerical simulations as well as formal asymptotic analysis 
we refer to [DGHll, HNV12]. 

As explained above, the underlying gradient structure ensures that the systems approaches - 
after a short initial transient regime with large dissipation - at time < to <^ 1 a state with small 
dissipation. Assuming both that remains small and that a changes regularly (i.e., on the time 
scale 1) for all times t > we can describe the dynamics for <C 1 as follows. 
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Figure 3: Moment weight for the definition of ^ with a* < u < 0. For a € (cr«, a*) and ^ -C 1, almost all 
of the total mass is concentrated in three narrow peaks located at X_(cr), Xo{a), and X+((j), but only the 
peaks at X± (a) are dynamically stable. For cr < a* and tr > a* , the mass is concentrated for ^ <C 1 in a single 
stable peak at X_{a) and X_^_{a), respectively. 

Formation of peaks The small dissipation assumption implies (see also Remark 19 below) that 
the moment 

m-= [ {H'{x)-a{t))''g{t,x)dx (4) 
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is also small, and we conclude that all of the mass of the system must be concentrated in narrow peaks 
located at the solutions to H'{x) = cr{t), see Figure 3. We can therefore (at least in a weak*-sense) 
approximate 

e{t, x) « Sx_{a{t)) for a{t) < cr*, Q{t, x) « Sx+{ait)) for o-{t) > a* (5) 

as well as 

0{t, x) ^ XI ^iii)^x_(a{t)){x) for a{t) e {a^, a*) , (6) 
ie{-,o,+} 

where the partial masses are defined by 

/X* rx» r+oo 

Q{t,x)dx, mo{t) := / 0{t,x)dx, m+(f) := / g{t,x)dx. (7) 
-OO J X* J Xt 

Notice that m_ (t) + mo (t) + m+ (t) = 1 holds by construction and that the moment ^ can be regarded 
as the formal limit of the dissipation as — )• 0. 

Thanks to (5), we have mo{t) « ~ for a{t) < a* and the dynamical constraint implies 

X_(cr(t)) Ri £{t), which determines the evolution of a. Similarly, with a{t) > a* we find m-{t) ^ 
mo{t) ^ and X+(cr(t)) £{t). These results reflect that H^^ is a single- well potential for both 
(T < cj* and a > a* attaining the global minimum at X_(cr) and X^{a), respectively. 

In the case of a{t) G {a^:, a*), the corresponding effective potential has two local minima and 
a local maximum corresponding to the three possible peak positions in (6). The peaks located 
at X±(^a{t)) are dynamically stable because adjacent characteristics of the transport operator in 
(FPi) approach each other exponentially fast for H"{x) > 0. Moreover, asymptotic analysis of the 
entropic effects reveals that each stable peak is basically a rescaled Gaussian with width of order 
v / \/ H" {X±{a)) . A peak at the center position Xo((t), however, is dynamically unstable because the 
spinodal characteristics separate exponentially fast with local rate proportional to r, and because the 
width of each peak is at least of order v. Each possible peak at Xo(t) therefore disappears rapidly, 
and by enlarging to if necessary we can assume that mo{t) ^ for all t >to. (This is different to the 
slow reaction regime, in which unstable peaks can survive for a long time due toO<i^<Cr<Cl). 

In summary, for any time t > to with a{t) G (o"*, a*) we expect that almost all of the mass is 
concentrated in the two stable peaks at X±{a{t)). In the limit — >^ 0, we therefore have mo{t) = 
and hence 

m-{t) + m+{t) =0, e{t) =m-{t)X-{a{t)) + m+{t)X+{a{t)) , 

where the last identity stems from the dynamical constraint. Notice that these formulas hold also 
for a{t) < cr* and a{t) > a* with m+(t) = and m-{t) = 0, respectively. 



Dynamics of peaks It remains to understand the dynamics of the multiplier a(t) and the partial 
masses m-{t) and m+{t) in the case of a{t) G (cr*, a*). The key observation is that although both 
peaks are spatially separated they can, at least in principle, exchange mass by a Kramers-type phase 
transition. In the many-particle picture this means that particles cross the energy barrier between 
the two wells of H^r due to stochastic fluctuations. Kramers investigated this large deviations problem 
in the context of chemical reactions in [Kra40] and derived his seminal formula for the effective mass 
flux between wells. In our notations, and with respect to our time scaling, this mass flux is, to leading 
order in v, given by 

-m-{t) « +m+{t) « m_(i)F_(t) - m+{t)F+{t), tF± « C±{a) exp (^-^^^^^ , (8) 

where the constants C±{a) do not dependent on v and are provided by Kramers' formula. In our 
context, however, the particular values of C±{a) are not important because the dominant effects 
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are the dynamical constraint and the time dependence of the energy barriers h± . For more details 
on Kramers' formula and the connection to the theory of large deviations we refer, for instance, to 
[HTB90, Berll]. 

We next discuss the implications of (8) for the different ranges of a. 

Subcritical case: For a{t) G (cr^, (t#) , the energy barrier between the two wells is larger than the 
critical value. This means 

h±{a{t)) > h#, F±{t) - exp « 

and we conclude that there is virtually no mass exchange between both peaks. In particular, the 
macroscopic dynamics reduces to 

m±{t) = 0, £{t) = &{t)(m-{t)X'_{a{t)) + m+{t)X'^{a{t) 

and describes that both peaks are transported by the dynamical constraint, see the right panel in 
Figure 4. 



X.{a{t)) 




X.{a{t)) 



Figure 4: Left panel: For supercritical a < a^, all the mass is contained in a single stable peak, which is 
located at X^{a) and transported by the dynamical constraint. (A similar statement holds for supercritical 
<T > cr^fi-) Right panel: For subcritical a € cr"^^, the mass is in general concentrated in two stable peaks, 
which are located at X- (a) and X+ (a) , and move according to the dynamics of i. Both panels: The width 
of each peak is proportional to u/^yH"{X), where X denotes the position, and the arrows indicate that the 
peaks can move to the left (for i' < 0) or to the right (for ^ > 0). The shaded regions in light and dark gray 
represent the intervals [X_(fT#), X+((t#)] and [x*, x*], respectively. 




X_(a#) 



X_((T#) 



X+{a#) 



Figure 5: For critical a, the coexisting stable peaks exchange mass by a Kramers-type phase transition, where 
C7 = (T# (left panel) and ct = a* (right panel) correspond to negative and positive mass flux, respectively. 



Critical cases: For a{t) = a^, we find 

h+{a{t)) = h#<h-{ait)), 



F_(t)«l~F+(t), 
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which means particle can move from the well at X-|-(o"(t)) towards the well at X_(cT(t)), but not 
the other way around. It is therefore possible that in the limit v ^ there exist time intervals of 
positive length with 

cr(t) = cr#, +m+{t) = -m-{t) < 0, £{t) = m^{t)X^{a#) + m+(i)X+ (c7#) . 

Similarly, it can also happen that the macroscopic dynamics is given by 

a{t) = a*, = -m_(i) > 0, t{t) = m-{t)X_[a*) + m+{t)X+{a*), 

reflecting an effective mass flux from the well at X_((7(t)) towards the well at X+((T(t)). Both 
critical cases are illustrated in Figure 5. 

Supercritical cases: For o'(t) G (c*, cr^), we verify 

h+{a{t)) < /i# < h-{a{t)), F-{t) < 1 < F+(t), 

and conclude that particles escape very rapidly from the well at Xj^{a{t)) but are trapped inside the 
other well at X-iait)). The only consistent choice for the macroscopic dynamics in this case is 

m-{t) = 1, m+(i) = 0, t{t) = X'_ {a{t))a{t) 

corresponding to the transport of a single stable peak, see the left panel from Figure 4. To be more 
precise, for states with a{t) G (a*, a^) and m^{t) > 0, the mass-dissipation estimates derived below 
imply that 'D{t) is large, and hence we expect that such states cannot be reached dynamically. (If 
such states are imposed in the initial data, a very rapid mass transfer during the initial transient 
regime produces m_(to) ~ 0.) Similarly, for a{t) G (a*, a*) the macroscopic evolution reads 

m-{t) = 0, m+{t) = 1, l{t) = X'^{a{t))&{t) 

and can be justified by analogous arguments. Notice that the limit dynamics in the supercritical 
cases is the same as in the single- well cases a{t) < cr* and a{t) > a*. 

2.4 Rate-independent model for the limit dynamics 

The above formulas for the limit dynamics can be translated into closed evolution equations for i, 
a, and the phase fraction n := m+ — ni-. These equations are illustrated in Figure 6 and are rate- 
independent because the macroscopic solution corresponding to i{t) = £{ct) with c > is given by 
a{t) = a{cf.) and jl(t) = fi{ct). For more details on the general theory of rate-independent systems and 
the different solution concepts we refer to [Miella]. Moreover, the limit dynamics exhibit hysteresis 
in the sense that the value of the output a at time t depends not only on the instantaneous value 
of the input £ but also on the history of the evolution (or, equivalently, on the state of the internal 
variable fi). 

In order to give a precise formulation of our limit model, we now define an appropriate notion of 
solutions. To this end we observe that the parameter constraints 

{X_(cj)} for £7 < a#, 

[X_(a), X+(a)] for ae a#] , (9) 
{X+{a)} for a>a* 



II G [-1, 1] , a G M, £e < 
define the macroscopic state space 



J7 := |(£, /i) G M-'^ satisfying (9)| (10) 

and that the dynamical constraint can be written as C{£, a, fj,) = 0, where the function 

Ci£, a, := ^X_{a) + ^X+{a) - £. (11) 

is well-defined and continuously differentiable on Q. We also recall that any Lipschitz function 
admits a classical derivative in almost all points (Rademacher's Theorem, see for instance [DiB02, 
Proposition 23.2] ). 
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/ 



■#.. 




a 



{m = +i} 



X* 



Figure 6: Cartoon of the macroscopic limit dynamics in the (£, (T)-plane. The gray soHd curve is the graph 

of iJ', the dashed black lines represent the level curves of /i, and the solid black lines correspond to the critical 
values (7^ and cr'^ , for which mass transfer according to a Kramers- type phase transition is feasible. The black 
and gray arrows indicate the evolution for decreasing and increasing £, respectively. Microscopic dynamics for 
small V. The evolution of q along the level sets of /U is illustrated in Figure 4, whereas the panels in Figure 5 
correspond to a{t) — and a(t) = cr*. 

Definition 10 (solutions to the limit model). A pair (a, ii) G C°'^([0, T] ;]R2) called a solution 
to the limit problem for given I G C°'^([0, T\), if the pointwise relations 



fi{t) = if a{t)i{a#,a*], < if a{t) = a#, > if a{t) = a* (13) 



hold for almost all t G [0, T] . 

In Appendix B, Proposition 34 we prove that for each i as in Assumption 5 and any admissible 
choice of the initial data (<t(0), //(O)) there exists a unique solution to the limit model, which is 
moreover piecewise continuously differentiable. We also emphasize that the limit model is equivalent 
to a constrained variational inequality. More precisely, introducing the convex functional 



Here, d means the set-valued derivative in the sense of subgradients, and the dynamical constraint 
enters via the pointwise relations (12). 

The heuristic derivation of the limit dynamics in §2.3 relies on two crucial assumptions for t > to, 
namely (a) that ^{t) is pointwise small, and (6) that &{t) is pointwise of order 1. In numerical sim- 
ulations one observes such a nice behavior but our convergence proof is based on weaker statements, 
which are, however, sufficient for passing to the limit — >■ 0. Specifically, below we only show (a) 
that ^(t) remains small, and (6) that a is Lipschitz continuous up to some small error terms. 

3 Auxiliary results 

The quantities c, C, and a always denote positive but generic constants (so their value may change 
from line to line) which are independent of but can depend on H, i, T, the constant from Assump- 
tion 6, and other parameters to be introduced below. Notice that the scaling law between r and ly, 
see Assumption 4, implies that a given positive quantity is exponentially small with respect to u if 
and only if it is bounded by Cr" for some constants a and C independent of v. 



{i{t), a{t), iJL{t)) e n with C{e{t), a{t), n{t)) = 
are satisfied for all t G [0, T], and if the dynamical relations 



(12) 




the dynamical relations (13) can be formulated as 

a{t) e d^n{^i{t)) + ^^x{^x{t)). 
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3.1 Mass-dissipation estimates 



In this section we derive mass-dissipation estimates, that means we show that small dissipation 
requires the total mass of the system to be concentrated near either both or one of the stable peak 
positions X-[a) and X^{a). These estimates become important in §4 because they guarantee (in 
combination with the L^-bound for V) that for each time ti there exists another time t2 € [ti, ti + t^^] 
with < /? < 1 at which the data are well-prepared. In the present section, however, all arguments 
and results hold pointwise in t and thus we omit the time dependence in all quantities. 

For the following considerations we introduce, for each cr G M, the relative equilibrium density 

7cr(x) := exp (^--^^^ > j '^"^^^ 

see Figure 7 for an illustration, and denote by 7cr,- and 7^ + the restriction of 7o- to the intervals 

I„- := (-00, Xoia)) and := {Xo{a), +00) , 

respectively. The functions 70- are naturally related to states with small dissipations. In fact, ja/za 
is the global equilibrium of the linear Fokker-Planck equation (FPi) with a{t) = a = const, and the 
energy functional 



£(t{q) '■= T^"^ I q{x)\u. q(x) dx + I H„{x)Q{x)dx 



just gives the relative entropy of g with respect to 7o-, that is 

£Aq) = 1^' [ Qi 

JR 



x) In ( 1 dx. 




iil/7a(a;) 



X_(a) Xo{a) X+{a) 




l/Eo 



X_((t) Xo{a) X+{a) 



Figure 7: The relative equilibrium density 7^^ for cr* < < 0, where Ej 



exp(-F,(X,(a))/i/2 



For 



1/ <C 1, the density 7o- exhibits two peaks located at X-{a) and X^{a). The width of these peaks is of order 
0{u) and the mass ratio between the peaks scales with exp ((/i_((t) — h+{a))/u^). The inverse density l/'ja 
forms a peak at Xo{a) and grows very rapidly for a; ^ ±00. 



3.1.1 On Poincare and Muckenhoupt constants 

We now summarize some well-known facts about L^-measures, which allow us to establish mass- 
dissipation estimates in §3.1.3. Within this section, let I = (a:_, be some (bounded or un- 
bounded) interval, 7 be a positive L^-function on the interval /, and Cp{'y) the Poincare constant of 
7, that means 

/ {iv'{x))'^j{x)dx / w{x)^{x)dx 
t — inf -h — il 

P\D weL {'ydx) / (yj(^x) — Wa,v^^^{x) dx / j{x)dx 
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where w' abbreviates the derivative of w with respect to x and L^(7da:) := {w : j^w{x)^^{x) d.x < 
00} . For each xq G /, we also introduce the one-sided Muckenhoupt constants C^(7, a;o) with respect 
to Xo by 



CS(7. X.) := ^^sup^^^ (/" ^ d,) (£ ,(y) iy) , 



It is known, see the discussion in [Fou05, Schl2], that 7 admits a finite Poincare constant if and only 
if the Muckenhoupt constants are bounded. 

Lemma 11 (Muckenhoupt constants bound Poincare constant). We have 

Cpi'j) < 4max|c^(7, xq), C^il, xq)^ 

for all 7 and any xq G /• 

Proof. The proof can be found in [Schl2, Proposition 5.21]. □ 
We also mention the lower bound 

Cp(7) > |max|c^(7, x^ed), C'm(7, x^ed)], 

where Xmed is the median of 7, which is defined by ff j(y) dy = ^ j(y) dy, and that Cf,, can 
easily be estimated for logarithmically concave functions 7. 

Remark 12 (C^ for logarithmically concave 7). Let 7(x) = exp ( — F(x)) . For convex and strictly 
increasing potential V : [xq, +00) — > M we have 

CiJ,(7,xo)<sup^. 



Similarly, the estimate 



ri- ( \ ^ X - Xo 
C^(7,xo)<sup^ 



holds provided that V : (—00, xo] — > M is convex and strictly decreasing. 
Proof. For each x > xq we estimate 

rx rx 

/ —(-\ dy = exp {V{y)) dy < exp (F(x)) (x - xo). 

Jxo Iky) J XQ 

Moreover, using Taylor-Expansion of F at x as well as the monotonicity of V' we find 

/•OO /"OO 

/ l{y)dy= I exp(-F(y))dy 

J X J X 

< exp ( - V{x)) exp ( - V'{x){y - x)) dy = ^""^ l^"^^""^^ • 

The first claim now follows immediately, and the arguments for the second one are similar. □ 

The mass-dissipation estimates derived below rely on asymptotic expressions for the Muckenhoupt 
constants of 7o- and the following observation. 
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Lemma 13 (variant of Poincare inequality). For any 7, the estimate 

/ w{xf^{x)dix < 2Cp(7) / {w'{x)) -f{x)dx + 2Cj{-f) / w{xf^{x)dx, Cj{-f) := ^ 

•'■^ •'^ ■'^ j 7(x) dx 

holds for all w G V?{^dx) and any subinterval J C. I. 

Proof. Thanks to 2ab < rja^ + ?7~^6^ and Holder's inequality we have 

and with rj := 2 Jj 7(0;) dx we find 

y (w(x) - Wav)^7(aj) dx > J (w{x) - i«av)^7(a^) d^; 

= / w{x)'^^{x) dx + w^y / 7(0;) da; — 2u;av / ty(x)7(x)da; 
J J J J J J 

> 2 / w(x)^7(a;) dx — / 7(0;) dx. 
J J J J 

Holders inequality also implies 

(x) dx 



J tu(x)^7(a 



x) dx 



and combining the latter two estimates with the Poincare estimate for w and 7 yields the desired 
result. □ 

3.1.2 Asymptotics of Poincare constants for 7^ 

In this section we derive upper bounds for the Poincare constants for 7cr and 7(t,±- For fixed a, the 
key observations concerning the z^-dependence can be summarized as follows. 

1. For C7 > C7* or C7 < (7* we find 

CpM = c^y'' 

because H^r is a regular single- well potential that grows quadratically at infinity. 

2. For cr G (<t*, ct*), the Poincare constant for small u is given by 

min|/i_(cr), /i_|_((j)}\ 



Cp(7(t) = C'^exp 



1/2 



because is a genuine double-well potential. This implies 

(a) Cp{'ja) ^ I/t for supercritical a as the energy barrier is smaller than the critical barrier 

= h-i^a"^) = /i+((T#), but 

(b) Cp{'ja) S> l/r for subcritical a since the energy barrier exceeds h^. 

Moreover, the Poincare constants of j^,- and 7^,+ are bounded by some constant C^- indepen- 
dent of u. 
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X_(a) 



Xo(a) X+(a) 



Figure 8: Examples of the effective potential H„ with cr < cr* (left panel) and cr* < cr < (right panel). 
Muckenhoupt constants ^^(70-, X-{a)) are estimated in the proofs of Lemma 15 and Lemma 16 



The 



For our purposes, however, we need estimates that hold uniformly in certain ranges of a and are 
derived in the subsequent series of lemmata. 

Lemma 14 (Poincare constants of 7(j,± if is a double- well potential). For each £ with < e < 
\{cF* — (7*) there exists a constant C, which depends only on e and H , such that 

Cp(7a,±) < C 
holds for all < 1/ < 1 and a € [a* + £, a* — s]- 

Proof. Let a e (a* + e, a* — s) be given. Since the potential H^r is strongly convex and strictly 
decreasing on the interval (—00, X-{a)], Remark 12 provides 

1 X — X_((t) 2 X — X_((j) 

< V sup — — = V sup 



< < 

mf a;<x_ (ct) -H"" (a^) mf a;<x_ {a* -s) H"{x) 
Moreover, Her is strictly increasing on the interval [X_((t), Xo{a)], and thus we estimate 

C^i(7.,-, ^-(^)) = sup ( r exp f +^') dy^ ( /"''"^'^ exp f-^") dy 



lup [ \ 

a),XoM] \JX 



E[X_(<7),XoW] 

sup {x-X-{a)){Xo{a)-x) <C. 

xe[X-(a),Xo(a)] 

Prom Lemma 11 we now conclude that 

Cp(7a,-)<max{Ci.2, C}, 
and the corresponding estimate for 70- + follows by symmetry. 



□ 



Lemma 15 (Poincare constant of ^ja- if H^r is a single- well potential). For each e > there exists a 
constant C, which depends only on e and H, such that 

Cp(7a) < C 

holds for all < u < 1 and a e [a* — 1/e, cr*] U [a*, a* -|- 1/e]. 
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Proof. Suppose that a G [cr* — cr*]. The potential is strongly convex and strictly decreasing 
on the interval (—00, X_((t)], and hence we show 



1/2 



C7^(7.,X_(a)) < — — - = Ci.' 



as in the proof of Lemma 14. In order to estimate 



we choose X > and notice that i^o- is strongly convex and strictly increasing on [X, +00), see 
Figure 8. In particular, we have 



x>X X — X x>X X — X 



with c := \rd^^j^H"{x) > 0, and Remark 12 yields 



x>X 

For a; > X we therefore obtain 



where 

' "JC 1 \ / /■+00 

Moreover, for x G [X_(cr), X] we estimate 

< {x-X-{a)){X -x)+Ca, 

where we used that 70- is strictly decreasing on [X- (a) , X] . Combining all estimates derived so far 
with Lemma 11 gives 

CpM<C{1 + u'')+C^, 

and thus it remains to bound C^. To this end we employ the monotonicity properties of H„ and H'^ 
to find 

The discussion in the case of cr G [a*, a* + l/e] is analogous. □ 
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Lemma 16 (Poincare constant of 70- if Ha is a supercritical double- well potential). For each e with 
< e < min{ci"# — cr*, a* — a*} there exist constants a and C which depend only on e and H such 
that 

holds for all a G [cr*, cr# — e] U [a* + £, cr*] and all sufficiently small u > 0. 

Proof. By symmetry and continuity, it is sufficient to consider the case a G (cr* , cr* — e] . This 
implies 

X_((t) < X* < Xo{a) < X* < X+((7) < X+((7#), 

and as in the proofs of Lemma 14 and Lemma 15 we derive the estimates 

,,2 



and 



where X := X+(cr#). Due to the monotonicity properties of and H'^, see Figure 8, we further 
obtain 

as well as 



'X_((t) 7^(2/) y \Jxo{a) 

< C3 exp 

We now abbreviate 



h+ia] 



and discuss four different cases: With a; > X we estimate 



For X e [X+{a), X] we find 

and since Her is strictly increasing on the interval [X+(cr), , there exists a constant C4 such that 



fa{x) <C4 1 + z^^ + exp 
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In the case of x G [Xo{a), X+(cr)] we verify 



and for x G [X_(cr), -'^o(c)] we finally get 



< C5 ( 1 + i/^ + exp 



where the last inequality holds since H^^ is strictly increasing on the interval [X_(cj), Xo(a)]. Taking 
the supremum over all x > X-{a) we now obtain, thanks to Lemma 11, the bound 



Mo) 

1.2 



Cp(7) < max{C- (7a,X_((7)), (7., < C (^1 + + exp 

The claim now follows with h+{a) < and since we have 

where o(l) means arbitrary small for small v. □ 
3.1.3 Estimates for the mass near the stable peak positions 

In order to establish the mass-dissipation estimates, we introduce the dissipation functional 

r (u''d^g{x) + {H'ix) - a)g{x)Y 
v,{e):= / ^ — ^-dx, 

and observe that 



T^(t{q) = ^^"^ I [dxvof'y^dx, / Qdx= / w^-^a^x for g = w^^a- 
Jr Jj Jj 



(14) 



Our first mass-dissipation estimate implies for each cr G (a* , a* ) that the mass is concentrated near 
the stable peak positions X_(cr) and -^^+(0-) provided that the dissipation is sufficiently small. 

Lemma 17 (upper bound for mass outside the stable peaks). For each e and any 77 with 

< £ < \{a* - (7*), < 7? < minjx* - X^{a* - e), X+(cr* + e) - x*| 
there exist constants a and C, which depend only on e and r], such that 

[ Qix)dx<CT'^('^^ + l) 

for all a G (a* + £, a* — s), any smooth probability measure g, and all sufficiently small v > 0. 
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Proof. Due to the bounds for ry, there exist constants C and a < 1 such that 



/ 7(j(a;)dx / 7^(0;) dx 

J I^,-\Br,{,X-{a)) ^ J I^,+\Br,{X+{a)) ^ 

J I<y.— J Ia.+ 



for all sufficiently small v. Using Lemma 13 twice with 

7 = 7o-,±, I = Ia,±, = qIIg-, J = Ia,±\Br,{X±{a)) 

we therefore arrive, see also (14), at the estimate 



f 



e{x) dx < 2(Cp(7.,-) + Cp(7.,+)) + Cr^ 



,{Br,(X_{a))UBr,{X+{a))) 

Moreover, the combination of Lemma 11 and Lemma 14 yields 

Cp{ja,±) < C, 

and this implies the desired result due to ^ w'^j dx < J^gdx = 1 and since we have i'~'^t < 
for all sufficiently small > 0. □ 

The second mass-dissipation estimate applies to strictly supercritical a and reveals that the 
dissipation controls the mass near the global minimizer of H„, which is given by X-{a) and X^{a) 
for a < and a > cr*, respectively. 

Lemma 18 (upper bound for mass outside the most stable peak). For each s > and any rj with 

0<V< {x* -X^{a* -e), X+ ((7# + e) - } 

there exist constants a and C, which depend only on e and rj, such that the implications 

a>a* + e =^ f gdx > 1 - Cr" + l) 

JBn{X+{a}) \ T J 

and 



a < (T# — £ 



[ gdx>l-Cr-(^ + l) 

JBniX.ia)) \ T J 



hold for any smooth probability measure g and all sufficiently small v. 

Proof. We only prove the first implication; the second one follows by analogous arguments. By 
Lemma 15 and Lemma 16, there exist positive constants C and a such that 

Making a smaller and C larger (if necessary) we can also assume that 



/ 

Jr\b 



-faix)dx 



■ya{x)dx 

for all sufficiently small v. The assertion now follows by applying Lemma 13 with 7 = 7^, / = R, 
and J = R\B^{X+{a)). □ 
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3.2 Dynamical stability of peeiks 

The most fundamental part of our convergence proof is to show that for sufficiently small v any 
solution to the nonlocal Fokker-Planck equation (FPi) + (FP2) can - at each sufficiently large time t 
and depending on the value of (j{t) - be approximated by either two or one stable peaks located at 
X_((T(t)) and/or X_|_((T(t)). In view of the mass-dissipation estimates from §3.1 it is clear that such 
an approximation is possible if the dissipation is small, but our approach lacks pointwisc estimates 
for Tyit) (we only have an L^-bound showing that Vit) becomes small after a small waiting time). 

We therefore control the approximation error by certain combinations of the moment ^ and the 
partial masses m_, mo, and m+, which all are defined in (4) and (7), because these quantities 
can be bounded pointwise in time. In order to identify the relevant combinations, we recall that 
m_ + mo + m+ = 1 holds by construction and that any solution to the nonlocal Fokker-Plank 
equation evolves according to the limit model if and only if 

1- C(0 + "^o(0 is small for all t, 

2. m+(t) for a{t) < a^, 

3. m+(t) Si 1 for a{t) > a*, 

4. m+(i) is almost constant for < a{t) < a*, 

5. m+(t) is essentially decreasing for a a^f, and essentially increasing for a ^ a'f^. 

In this section we derive upper bounds for ^ (t) + mo (t) and discuss the evolution of m_ and m+ 
afterwards in §3.3 and §3.4. We start with some auxiliary results which hold pointwise in time and 

does not rely on dynamical arguments. 

Remark 19 (dissipation bounds There exists a constant C such that ^{t) < T>{t) + Cv^ for all 

t>0 and v > 0. 

Proof. The definition of V, see (2), implies 

Vit)= [ ({H\x)-afg{t, x) + z/^ ^^"f/' "^^^ + 2i.^{H'ix) - a)d,Qit, x)] dx 
JR \ Q\t, X) j 

> [ {H'{x)-afQ{t,x)dx-2u'^f H"(x)Qit, x) dx, 
Jr Jr 

and this gives the desired result since we have \H"{x)\ < C for all a; G M. □ 

Lemma 20 (error terms in algebraic relations between £, a, and m±). For each e with < s < 
|(cr* — (T*) there exists a constant C , which depends on e but not on v, such that the implications 

|£(t) - X_ {ait)) I < C^y^{t)+mo{t)+m+{t) 

as well as 

ait) G[a,+e, a* -e] =^ |£(t) - m_(i)X_((7(t)) - m+(i)X+((7(t))| < C ^it) + moit) 
and 

ait) G [a* + £, +oo) =^ |^(t) - X+ {ait)) \ < Cv'C(t) + m_(t) + mo(t) 
holds for all < t < T and all v > 0. 
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Proof. We only prove the first impHcation; the derivations of the second and the third one are similar. 
By definition of the partial masses rrii, i G {—,0,+}, we find 



/X r+oo 
{x-X_{ait)))Q{t,x)dx+ XQ{t,x) 
-OO J X* 



+00 

dx 



and hence 



(r-x* r+oo \ 

mo{t) + m+{t) + J \H'{x) -a{t)\g{t, x)dx + J XQ{t,x)dxL 

where we used \a{t)\ < C, the asymptotic grow of H', and that |a; — X-{a{t))\ can be bounded by 
\H'{x) — a{t)\ because of 

0<c<H"{X-{a)) <C <oo for all a e {-00, a* - e] . 

Moreover, Holder's inequality yields 

rx \ ' 

m-{t) I \H'{x)-a{t)\^Q{t,x)dx\ < 
J —00 J 

thanks to m^{t) < 1, as well as 

r+oo 1/2/ r~^^ \ 1/2 

/ xg{t, x)dx < (rno{t) + m+{t)) ( / x^g{t,x)dx\ < C\/mQ{t) + m+(t), 



thanks to f^x'^g^t, x)dx < C. The first implication now follows from combining all result (notice 
that mi{t) < y/mi{t)). □ 

The assertions and the proof of Lemma 20 can easily be generalized to other moments. 

Remark 21. For any continuous moment weight ip that grows at most linearly and each e as in 
Lemma 20 there exists a constant C, which depends on e and tp but not on v, such that 



i;ix)g{t,x)dx- mj{t)i;(^Xj{a{t))') 

je{-,+} 



<C^m + moit) 



holds for all sufficiently small v as long as a{t) G [c* + e, a* — e]. Moreover, similar results hold in 
the cases a{t) G (—00, a* — e] and a{t) G [a* + e, +00). 

3.2.1 Evolution of the moment ^ 

We next study the dynamics of ^ and derive an upper bound for ^(t) by combining the moment 
balance for ^ with a simple ODE argument. For the formulation of this result we define 



fx^+ri 
I x*—r] 

for all 77 > 0. 



rx,+r) 

mj^{t) := I g{t,x)dx 

J x*—r] 



Lemma 22 (pointwise estimate for For each ij > there exists a constant C , which depends on 
77 hut not on v, such that 

sup i{t) < + C( + sup mr,{t) ] 
te[ti,t2] \ te[ti,t2] J 

holds for all <ti < t2 <T and all sufficiently small v > 0. 
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Proof. Using the abbreviation tjj{t, x) := {H'{x) — CF{t))^ as well as (FPg) and integration by parts, 
we easily compute 

T^{t) = -2T&{t) / {H'{x) - a{t))g{t, x)dx+ / tpit, x)Tdtg{t, x) dx 

= +2T^&{t)eit) + [ x)d^(i^^d^Q{t, x) + {H'{x) - a{t))Q{t, x)) dx 
Jr. ^ ' 

= +2T'^&(t)i{t) + 1^^ ip"it, x)git, x)dx-2 / H"{x)ip{t, x)Q{t, x) dx, 

JR JR 

as well as 

Ta{t) = riit) + f H"'{x)git, x)dx- [ H"{x){H'{x) - a{t))g{t, x)dx. 

JR JR 

In view of 

x)\ + \H'{x)\ + \H"{x)\ + \H"'{x)\ < C(l + x^) 

and 

+ K'WI + Wit)\ +/"(! + x^)Qit, x)dx<C 

JR 

see Assumption 1, Assumption 5 and Lemma 7, we therefore find 

riit) < C{y^ + r) - 2 j H"{x)'4j{t, x)g{t, x) dx. 

JR 

Moreover, since H is smooth, there exist constants c and C such that 

H"{x)>c for all x e R\[x* - i], x^ + rj], \H"{x)\<C for all x e [x* - r], x^ + -q] , 
and this implies 



/ H"{x)i/j{t, x)g{t, x)dx> / H"{x)il;{t, x)g{t, x)dx+ / H"{x)il;{t, x)g{t, x)dx 

JR J J 

R\[x*-r),x,+ri] x*-v 

>c J ipit, x)g{t, x)dx - Cnir^it) = c^{t) - CMr,, 



Xx*-r],x»+T]] 

where is shorthand for sup^g^^^ j^] ^??(^)- Consequently, find 

rm<-cm + C{iy^ + r + Mr,) 
for all t e [ti, t2], and Gronwall's Lemma finishes the proof. □ 

3.2.2 Conditional stability estimates 

We are now able to investigate the dynamical stability of peaks. More precisely, assuming that 
a{t) remains confined to certain intervals we now derive estimates that control the evolution of 
^{t) + mo{t). In the proof we employ, apart from the estimates for ^, local comparison principles for 
linear Fokker-Planck equations in order to show that only a very small amount of mass can flow into 
the unstable interval {x*, x*). 
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Xo{a) 



X+{^) 



Figure 9: Cartoon of the supersolution g and the characteristic lengths from the proof of Lemma 23; the 
hatched regions indicate intervals of length 2r/. The strictly decreasing and increasing branches are given by 

rescaled equilibrium solutions corresponding to If and tj, respectively. For < ^ 1, ^ is therefore very small 
in [x* — 77, a;* + rj] and exhibits boundary layers with width of order v near X-{a) and X^{a). 



Lemma 23 (first conditional estimate for ^ + mo). For each e with < e < ^{a* — a^) there exists 
a positive constant C, which depends only on £ but not on u, such that the implication 



a{t) G [cr* + £, cr* — e] for all t E [ti, t2] 



sup {^{t) + mo{t)) < cU{ti) + mo(ti) + 1^^) 
telti,t2] ^ ' 



holds for all t* <ti <t2 <T and all sufficiently small v > 0. 

Proof. Within this proof we regard (FPi) as a non-autonomous but linear PDE for g, that means 
we ignore (FP2) and regard a as a given function of time. 
Step 1: We first choose a, a G (a*, a*) such that 

= h-{a) = I min |/i+((7* + e), h-{a* — e)}, 

and the monotonicity properties of h- and /t+, see Figure 2, ensure that a G (a*, a* + e) and 
a G {a* — £, a*). Employing the monotonicity of X-, Xq, and X+, we easily verify, see also Figure 
9, the order relation 

< X_{a) <x*< Xoicr) < Xo(a(ti)) < Xo{a) < x, < X+{a) < X+{a{ti)) 

and thus we can choose rj > such that the distance between any two adjacent points in this chain 
is larger than 27/. In particular, by definition of ^ we find 



g{ti,x) dx<C(e(ti)+mo(ti)) 



(15) 



for some constant C depending on ij. 

Step 2: We define a local supersolution 'q on the interval [X_(a), X_|_(a)] by combining rescaled 
versions of the monotone branches of 7^ and 7^. More precisely, we set 



q{x) := V ^ < 




^ —] for X-{a) < X < Xo{a), 



H^{X+{a)) - H^{x) 



for Xo(ct) < a; < Xo(o:), 
for Xo(a) <x<X+{gi). 



Our choice of a and a implies that g is continuous with 

-g{X_{a)) =g{X^{a))=y-\ 
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and thanks to our choice of 77 we readily show that there exist constants a and C such that 

'q{x) < Ct"' for all x e [x* — r), + 77] 

and all sufficiently small u. Moreover, g is by construction continuously differentiable and piecewise 
twice continuously differentiable, where 



d^(u^d^Q{x) + {H'{x) - a{t))Q{x)) 



' (a - a{t))d^:Q{x) for <x< Xo(a), 

H"{x)q{x) for Xo(a) < ar < Xo((t), 

. (a - a{t))d:,Q{x) for Xq{o) <x< 



Combining this with 

a<G{t)<G for iG[ii,t2], < for x G [Xo(a), Xo(o:)] 

and 

d^Q{x)<Q for xG Xo(a)], ^^^(a;) > for x G [Xo(a), 

gives 

+ {H'{x) - < = TdtQ{x), 
and we conclude that ^ is in fact a weak supersolution to (FPi) on the space-time domain [ti, i-2\ x 

S'tep 5; We consider three solutions qq^ and ^+ to (FPi) on the time interval [ii, defined 
by the initial conditions 

x) = x)x(-oo,x_(a)](a;), 
£»o(ti, a;) = x)x(x_(a),x+{a)](2;), 
x) = x)x[x+(o:),+oo)(a;), 

where x/ is the usual indicator function of the interval /. All three functions are nonnegative, and 
thus we find 

C 

Q±{t, X) < Q{t, X) <^ 

for all X G M and t>ti>U thanks to the L°°-estimates from Lemma 7. We now conclude that 

g±{t, X^{a)) < Q{X^{a)), e±{t, X+{a)) < 

for all t £ [ti, t2], and the comparison principle yields g±{t, x) < g{t, x) for all t E [ti, ^2] and almost 
ah X G [X-{a), X+{a)]. We therefore get 

/ {g-{t, x) + g+{t, x)) dx < 2 / ^(x) dx < Cr" 

Jx*—r] Jx*—r^ 

for all t G [^1, ^2]- On the other hand, using the mass conservation property of (FPi) we estimate 

rx*+r] r+00 r^+(z) 

/ ^o(^)2;)dx< / go{ti,x)dx= / ^(ii,x)dx. 

Jx*-r] J -00 Jx-{a) 

With g = g- + go + g+, the estimate (15), and by taking the supremum over t we therefore get 

r,{t) <c(^J^ ^^^^ g{ti, x) dx + r'*^ < C(^{ti) + mo(ii) + u^) , 
where we used r" < i^^. Finally, the desired result follows from Lemma 22. □ 
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sup rrirj 
te[ti,t2] 



Lemma 24 (second and third conditional estimate for ^ + mo). For each e with e > there exists 
a positive constant C , which depends only on e hut not on v, such that the implications 

(T{t) > (7* + £ for all t G [ti, t2] =^ sup (^(i) + mo(i)) < C{^{ti) + mo(ti) + m_(ti) + u'^) 
and 

ait) <a*-e for all t G [ti, t2] =^ sup {^{t) + mo{t)) < cU{ti) + mo(ti) + m+{ti) + z/^) 

te[ti,t2] ^ ' 

hold for all t* < ti < t2 < T and all sufficiently small u > 0. 

Proof. The proof is very similar to that of Lemma 23, and thus we only sketch the main ideas. For 
the first implication, we set 



:= (T* — ^s, V '■— 2 



\ min ja;* - Xo(a:), - x*, X+((T* + e) - 



which in turn implies 

rX+{.a) rX+ia) 

/ g{ti, x) dx = m-{ti) + mo{ti) + ^(ti, x) dx < m_(ii) + mo(ii) + C^(ii) 

J —oo J 

for some constant C, which depends on rf and hence on e. We then define a local supersolution g in 
the interval [x*, X^{a)] by 

^ I exp(-^^±^J for x*<x<Xo{a), 

exp I = — 1 for Xo{a) <x< X+[a), 

and consider two solutions and g+ to (FPi) with 

g-oih, x) = g{ti, x)x(-oo,x+{a)]{x), Q+ih, x) = g{ti, x)x[x+{a),+oo){x)- 

Employing the comparison principle with respect to the space-time domain [ti, t2\ x [x* , X^{a)], we 
then show that 

nT'Tjit) = / g-o{t,x)dx+ / g+{t,x)dx 

Jx*—ri Jx*—ri 

/+00 pxt+rj 
g-Q{ti,x)dx+ / 'g{x)dx 
-oo Jx*—r) 
rX+{a) 

< / g{ti, x)dx + CT" 

J —oo 

< C(e(ii) + + mo(ti) + t"^ , 

and the assertion follows from Lemma 22. The second implication can be proven analogously. □ 

3.3 Mciss transfer between the stable regions 

We next investigate the evolution of m_ and m+ by means of appropriate moment balances. The 
resulting estimates imply for <C 1 that the mass flux from the left stable interval (— oo, x*] towards 
the right one [x*, +oo) is - up to small correction terms - positive for a{t) > but negative for 
a{t) < C7*, and hence that there is essentially no mass transfer in the subcritical regime a{t) G 
(cr#, cr'^). These findings perfectly agree with the large deviations results that we obtained in §2.3 
by analyzing the orders of magnitude for the different terms in Kramers formula (8) . 
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Figure 10: Cartoon of the moment weight V that is used in the proof of Lemma 25. The strictly decreasing 

branch of i]) on [x*, X_|-((7)] is given by the rescaled and shifted primitive of —1/70- and has effective width of 
order v. For ^ 1, the function i/' is therefore close to the indicator function of (—00, X^ia)^. 



Lemma 25 (almost-monotonicity estimates for m±). For each e with 

< £ < min{(7* — (T#, £7* — (T*} 

there exist constants a and C, which only depend on e, such that the implications 



a{t) > cr^ + £ for all t e [ti, t2] 



and 



u{t) <u* -£ for all t G [ti, t2] 



sup m_(t) < m_(fi) + mo(ti) + Cr" 



sup < m+(ti) + mo(ti) + Cr° 

ie[ii,t2] 



hold for t* < ti < t2 < T and all sufficiently small v > 0. 

Proof. We demonstrate the first implication only; the second one follows analogously. In what follows 
we control the evolution of an upper bound for m_ , namely the moment 



m. 



.(t) := / tlj{x)Q(t, x) dx. 



Here, the weight ip is defined as piecewise constant continuation of an appropriately rescaled and 
shifted primitive of —1/7^, where a is shorthand for a# + e. More precisely, we set 



i){x) := < 



1 



exp 



fx* 




exp 



dy 



for X < X*, 



for X* < X < X+{a), 



for X > X+{a), 



and refer to Figure 10 for an illustration. In particular, ip is continuous as well as piecewise twice 
continuously differentiable, and thus we readily verify (using (FPi) and integration by parts) the 
moment balance 



Tm-{t) 



i^'{x)(u^d^g{t, x) + {H'{x) - a{t))g{t, x)) dx 
= - i^'{x)[u^d^g{t, x) + {H'{x) - a{t))g{t, x)j dx 

■^"{x) + {a{t) - H'{x))ilj'{x)jg{t, x)dx. 



b.t. + 
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Here, the boundary terms are given by 

b.t. = u'^^'{x* +0)g{t, X*) - u'^^'{X+(a) - 0)^(t, X+{a)), 

and the notation ±0 indicates that the boundary values must be taken with respect to the interval 
[x*, X+{a)]. 

It remains to estimate m-{t) for all t e [ti, t2], that means for a{t) > a. We first infer from the 
definition of that 

il;'{x) < 0, u^2l;"{x) + {a- H'{x))il;'{x) = for ah x G [x*, X+{a)] . 

We next observe that due to h+{a) < hj^{(j^) the asymptotic properties of 7^ imply 

sup '^{x) < Ct'^, u-^ \tP'{X+{a) - 0)1 < Cr^+" 

X>Xt, 

for some positive constants a, C and all sufficiently small v. Finally, we have Q{t, X+(a;)) < Cu~^ 
according to Lemma 7. Combining all these estimates wc therefore find 

in-{t) < u'^\^'{X+{a) - 0)\Q{t, X+{a)) < Ct" 

for all sufficiently small u, and hence 

sup m-{t) < sup m-{t) <rn-{ti) + Ct'' <m-{ti) + rno{ti) + Ct°', 
te[ti,t2] te[ti,t2] 

where we used that t2 <T and m_(ti) < m_(ti) + mo(ti) + sup^^^^ '^(x). □ 
3.4 Monotonicity relations between a and £ 

As an important consequence of the almost-monotonicity relations for m_ and m+ we now establish, 
up to some (small) error terms, monotonicity relations between i and a. These results have three 
important implications, which can informally be summarized as follows: 

1. If a{t) ^ or a{t) holds for all t in some interval [ti, t2], then £ must be essentially 
decreasing or increasing, respectively, on this interval. The dynamical constraint then implies 
in the limit — )■ that the phase fraction ji is decreasing and increasing for a = and 
a = cr*, respectively. 

2. If [ti, is some time interval such that a behaves nicely with 

(a) (crossing from above) a{t2) < cr^ < cr{ti) < cj*, or 

(b) (crossing a* from below) < o"(ti) < a* < a{t2), 

then t2 — ti can be bounded from below by \o'{t2) — a{ti)\. This implies, roughly speaking, that 
solutions for small v cannot change too rapidly from subcritical a to supercritical a. 

3. If [ti, 12] is some time interval such that a stays inside the subcritical range (c#, c*), then 
\cr{t2) — a{ti)\ can be bounded from above by 1-^(^2) — £{ti)\, and this gives rise to Lipschitz 
estimates for subcritical a in the limit z/ — t- 0. 

Lemma 26 (conditional monotonicity relations). Let e be fixed with 

< £ < ^ minjo"* — o"#, cr* — o"*}. 

Then the implications 

a{t) e [a^ + e, a* -e] for all t e [ii, ^2] =^ 9{cT{ti) - a{t2)) < i{ti) - i{t2) + e.t. 
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and 

a{t)e[u^ + e,u* -e] for all te[ti,h] =^ g{cr{t2) - a{ti)) < £{t2) - i{ti) + e.t. 
hold with error terms 

e.t. := c(^^/^{tl)+mo{tl) + ^^(^2) + mo(t2)) + Ct" 

for all t-f < ti < t2 < T and all sufficiently small u > {). Here, g is the increasing and piecewise 
linear function g{s) = Csgii(s)S, where C_, C+ > are independent of both s and v, and a, C denote 
two constants which depend on e hut not on v. 

Proof. We only derive the first implication; the arguments for the second one are similar. For the 
proof we set x±{t) := X±(^a{t)) as well as 

i{t) := m^{t)x-{t) + m+{t)x+{t), e := VC(ii) + ^0(^1) + \/C(^2) + mo(t2), 

and suppose that a{t) G [cr* + e, a* — e\ holds for all t G \ti, t2\. Lemma 20 yields \i{ti) — l{ti)\ < e, 
and we conclude that 

Ce + i{ti) - l{t2) > i{ti) - i{t2) 

= E K(ii)-m,(t2))x,(ti)+ rnj{t2){xj{ti)-Xj{t2)). (1^) 

Thanks to m-{t) + mo{t) + m+{t) = 1 we find 

m_(ti) - m-{t2) = -(mo(ti) - mo(t2)) - {m-+{ti) - m+{t2)) 

and hence 

{mj{ti) - mj{t2))xj{ti) > (m+(ti) - m+{t2)) {x+{ti) - x-{ti)) - Ce, 

i6{-,+} 

where wc used that |x±(ii)| < C and mo(ti), mo(t2) < e. Moreover, Lemma 25 provides constants a 
and C such that 

m+{ti) - m+{t2) > -C{e + r") 
holds for all sufficiently small u, and in view of x+{ti) > x-{ti), see Remark 2, we arrive at 

J2 K(ti)-m,(t2))x,(ii)>-C(e + r"). (17) 

On the other hand, we have Xj(ti) — Xj{t2) = X'-{dj)(a{ti) — o{t2)) for some intermediate value dj, 
and the monotonicity properties of X_ and (again Remark 2) ensure the validity of 

Xj{ti) - Xj{t2) > g{cr{ti) - cr{t2)), 

where 

C_ := max max X'Aa), C4. := min min X'Ja), I := \h(a^ + ajt), J((7* + cr*)l. 

^e^ie{-,+} ^ ^e/ie{-,+} 

We therefore find 

^'jit2){xj{ti) - Xj{t2)) > g{a{ti) - a{t2)) (m_(ti) + m+(ti)) 

ie{-,+} 

> g{a{ti) - a{t2)) - Ce, 

and the desired implication follows by combining this estimate with (16) and (17). □ 
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4 Passage to the limit p ^ 



The arguments used in §3 to characterize the dynamics of the partial masses m_ and m+ depend 
cruciahy on the range of a. In particular, we have have quite strong results for the subcritical case 
a € o'*) because here both m_ and m+ are, to leading order in u, constant in time. Wc can 

also control the evolution in the supercritical cases a G (—00, a^) a G (a^, +00) because then cither 
m_ or m+ is always very small. In the critical cases a ^ cr^ and a ^ cr*, however, we can use only 
relatively weak monotonicity relations, and this complicates the analysis of the limit 1/ — >■ 0. Our 
strategy is therefore as follows. We introduce a small parameter e with 

< £ < := 2 min{cr* — cr*, cr# — cr*, a* — (7#} 

and accept to have only incomplete control over the dynamics as long as a{t) is inside the £- 
neighborhood of either cr^ or o"*. Afterwards we pass to the limit u ^ along sequences (en)„gN 
and (fn)neN with e„ — )• and < Vn < i^i^n) — ^ 0, where the critical value ^'(en) will be identified 
below. 



4.1 Approximation by stable peaiks 

Heuristically it is clear that the small parameter dynamics evolves according to the rate-independent 
limit model if and only if the state of the system can be approximated 

1. by two narrow peaks located at and X^{t) as long as a{t) G ((T#, <j*), 

2. by a single narrow peak located at X-{t) or X^{t) for a{t) G (—00, a^) or a{t) G (cr#, +00), 
respectively. 

In this section we combine all partial results from §3 to show that these assertions are satisfied for 
all sufficiently small u and e, and all sufficiently large times t. Specifically, we consider 

m+(t) for a(t) G (—00, — e] , 
for a{t) G (cr# — £, cr* + e) , 

m-{t) for a{t) G [cr* + £, +00) , 

and prove that is small for all times t >ti provided that the dissipation is small at time ti. This 
conclusion is in fact at the very core or our approach as it allows us to convert the L^-bound for the 
dissipation into moment estimates that hold pointwise in time. 



at) := at) + moit) + < 




Figure 11: Schematic representation of the intervals Ji as well as the cr-domains for the cases Nj and Pj as 
used in the proof of Lemma 27. For £ ^ 0, we have J_2 {—00, cr#], Jo — >■ [cr#, cr*] , J+2 — >■ [c*, +00) as 
well as J_i and J+i — )• {ct*}. 



Lemma 27 (pointwise upper bound for Q. For each e G [0, £*] there exist 'positive constants /3 < 1 
and C, which depend on s but not on v, such that the implication 

V{to) <T^ =^ sup C{t) < Cu'^ (18) 
te[to,r] 

holds for all t* <to <T and all sufficiently small u > 0. 
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Proof. We consider the intervals 

J_2 = (-00, cr# - e] , J_i = [(7# - e, (7# - ie] , 
as well as Jq := [a^ — a"^ + ^s] and 

J+i = a* + ^e, a* + e , J+2 = o-* + £, +00^ ■ 

These intervals and the different cases considered within this proof are illustrated in Figure 11. 

Part 1: We first prove the assertion under the assumption that a remains confined to at most 
two or three neighboring intervals from {J_2, J-i, Jq, J+i, J+2}, and start with the case 

cr(t) G J_2 U J_i for all te[ti,t2], {N_) 

where tQ < ti < t2 < T. Under this assumption, Lemma 24 combined with Lemma 25 provides 

constants ai and Ci such that 

sup {^{t) + mo(i)) < Ci Uiti) + + mo(ti) + 

as well as 

sup m+(t) < m+(ti) + mo(ti) + Cir"^ 

te[ti,t2] 

Moreover, by Lemma 18, there exist constants a2 and C2 such that 

mo(ti) + m+(ti) < C2r"^ (r-'P(ii) + l) , 
and Remark 19 yields a constant C3 such that 

We now choose Pi G (0, 1) sufhciently large such that 02 + /3i — 1 > and this guarantees (via 
^(ii) < (C'a + and m+(ti) + mo(ti) < u^), that the implication 

(iV_) and V{ti)<T^^ =^ sup (^(t) + mo(i) + m+(t)) < Qz^^ (19) 

holds for all sufficiently small > 0, where C4 := Ci (C3 + 3) . 
The arguments for the case 

a{t) G J+i U J+2 for all t G [ti, 12] (iV+) 

are entirely similar. In particular, possibly changing all constants introduced so far, we readily 
demonstrate that 

(iV+) and P(ii)<T^i =^ sup {i{t) + m-{t) + mo{t)) < dv"^ (20) 

holds for all sufficiently small > 0. 
We next study the case 

a{t) G J_i U Jo U J+i for aU i G [ii, t2] , {Nq) 

and observe that Lemma 23 provides a constant C5 such that 

sup {^{t) + mo{t)) < C5 i^ih) + moih) + u"^) . 
te{ti,t2) 
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By Lemma 17 we find further constants ae and Ce such that 

mo(ti) <C6T°«(r-ip(ti) + l), 

and we choose [52 G (0, 1) sufficiently close to 1 such that ag + /32 ~ 1 > 0. This ensures (via 
< (C's + l)i^^ and mo(ti) < z^^) that the implication 

{No) and X)(ti)<r'^2 ^ sup (^(t) + mo(t)) < Cyi^^ (21) 

«e[ti,t2] 

holds for all sufficiently small z/ > 0, where C^ := C5(C3 + 3). 
Part 2: We set 

/3 := max ^2}, C := max {C4, C7}. 

Our next goal is to demonstrate that whenever the systems passes for t G [ts, C [f*, T] through 
the entire interval J±i, then there exist at least one time t in between tj, and t4, at which the data 
are well prepared in the sense of P(i) <t^ . To this end, we have to discuss the four cases 

a{t) G J_i for all t G [ts, t/^ , a{t3) = a#-e, a{t4) = a^-\e (P_o) 

and 

a{t) G J-i for all t G [t^, U] , (7(^3) = (t# - \e, a{U) = - e (Pq-) 

as well as 

a{t) G J+i for all t G [ia, , (^(^s) = a* + £, a(i4) = cr* + (P+o) 

and 

(7(t) G J+i for all t G [t-^, U] , ^(ta) = a* + ^e, cj(t4) = a* + e (Po+) 

but by symmetry it is sufficient to study (P-o) and (Pq-)- We first discuss the case (Pq-) and 
suppose that 

C{t3) = ^it3) + mo{h)<Ci^^ 
Lemma 26 combined with the uniform bounds for \i{t)\ yields constants cg and Cg such that 

^4 - ^3 > C8(a(ia) - o-(i4)) - Cqv = ^cse - Cgv, 



and hence there exists a positive constant cio such that t4 — ta ^ cio for all sufficiently small v. Since 
we have j^^ P(t) dt < Cut, there exists at least one time t G [ta, t4] (which depends on u) such that 

V{t) <^T<T0 
ClO 

for all sufficiently small v > Q. We have thus proven that the implications 

(Pot) and C{h) < Cu'^ =^ V{t) < for some t G [ta, ^4] (22) 

hold for all sufficiently small > 0. 
For the case (Pq-) we assume that 
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Similar to the above discussion of the case (iV_), we use of Lemma 24 and Lemma 25 to show that 
there is a constant C12 such that 

CiU) = ^{U) + moiU) + m+iU) < Ci2Z/^ 

holds for all sufficiently small v > 0. Prom Lemma 20 we further infer that there is a constant C13 
such that 

\X_{a{t4)) - X_{a{t3)) \ < \£{U)-i{t3)\ + Ci3iy\ 

and the properties of X_ and i imply that t4 — is > C14 holds for all sufficiently small u > and 
some constant C14 independent of e. In particular, using J^^ dt < Cut once more, we show that 
implications 

(Pzpo) and C(*3) < Ciy"^ =^ V{t) < for some t G [^3, U\ (23) 

holds for all sufficiently small v > Q. Finally, we recall that we have 

(Pot) or (P^o) =^ U-t3>c (24) 

for some constant c > and all sufficiently small z> > 0. 

Part 3: Wc finally return to discussing the time interval [to, T] and show in a preparatory step 
that there exists a sufficiently large time Iq G (to, T) such that 

sup C(*) < Cv'^, V{to) < T^. (25) 
te[to,to] 

Suppose at first that cr(to) € J-2- If c(i) € J_2 U J_i holds for all t G [to, T\, then we are done with 
to = as (19) implies (18). Otherwise we consider the times 

t4 := inf |t G [to, T] : a{t) = cr# — ^e|, t^ := sup |t G [to, t4] : a{t) = cr# — e|, 

which are well-defined as a is continuous. By construction, the intervals [to, ts] and [ts, t4] corre- 
sponds to the cases (iV_) and (P-o), respectively, and the existence of to G [t3, t4] is a consequence 
of (19) and (22). Similarly, the case cr(to) G J+2 can be traced back to the cases (N^) and (P+o), 
and to is provided by (20) and (22). 

Now suppose that cr{to) G Jo. If <T(to) G J_i U Jo U J+i holds for all t G [to, T], we set to = T and 
are done by (21). Otherwise we find times ts < t4 such that [to, t3] corresponds to (Nq) and [t3, t4] 
to either (Po-) or (Po+), and the existence of to is implied by (21) and (23). 

For cj(to) G J±i, wc are either done via a{t) G J±i for all [to, T], or we find a time ti with 
a(t) G J±i for all t G [to, ti] and either a{ti) = ± e or cr(ti) = it ^e. Depending on the value 
of cr(ti), we can now argue as for a{to) G J±2 or c7(to) G Jo. 

In summary, we have now proven the existence of to with (25). Our arguments can easily be 
iterated, and since (24) provides a lower bound for the time covered by two subsequent iterations, 
we finally arrive at (18). □ 

4.2 Continuity estimates for a 

As a further key ingredient to the derivation of the limit dynamics we now show that a is, up to 
some error terms, globally Lipschitz continuous in time. These estimates become important when 
establishing the limit — t- because they imply the existence of convergent subsequences as well as 
the Lipschitz continuity of any limit function. 

Lemma 28 (Lipschitz continuity of a up to small error terms). For each e G [0, e*] there exist 
constants a and C , which depend on e but not on u, as well as a constant Co, which is independent 
of both e and u, such that 

|a(t2)-C7(ti)| <Co(|t2-ti|+e) +C(t"+ sup vW) 

te[ti,t2] 

holds for all t* <t\ <t2 <T and all sufficiently small v. 
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Proof. Step 0: We introduce appropriate cut offs in cr-space. More precisely, we define 
as well as 

where the nonlinear projectors P{a,a) given by P(^^ :^-^(a) := max{min{(7, ct}, a}. These definitions 
imply 

+2 

J2 ajit) = ait) + 2{a*-a#), (26) 

and since a is (for any given > 0) continuous in time, all projected functions aj depend continuously 
on t as well. 

Step 1: To show that ctq is almost Lipschitz continuous, we assume without loss of generality 
that ao{ti) < ao{t2) and consider at first the special case of ao(t) = a(t) G [(To(ii), 0-0(^2)] for all 
t e [ti, t2]- Under this assumption. Lemma 26 provides constants a, C and Co such that 

ko(i2)-ao(ii)| <CoK(t2)-^(ii)| + c(vW+vW + T") (27) 

holds for all sufficiently small u > 0. In the general case, we introduce two times ti and t2, which 
both depend on v, by 

ii := max {t G [ti, t2] : ao{t) = ao{ti)}, h := min {t G [ii, t2] : aoit) = ao{t2)}, (28) 

and notice that the Intermediate Value Theorem (applied to the continuous function ctq) ensures that 
(To is a bijective map between the intervals [ti, t2] and [ao{ti), cro(t2)]- In particular, our result for 
the special case applied to the interval [ii, £2] combined with \t2 — ti\ < 1^2 — ^i| yields again (27). 

Step 2: We next derive a Lipschitz estimate for o"+2- As above, we suppose that (7+2(^1) < c+2(i2) 
and consider at first the special case of (y+2{t) = a{t) G [17+2(^1), C7_|_2(i2)] for all t G [ii, t2\. From 
Lemma 20 we then infer that 

|%)-X+(a+2(ii))| <CvW, i = l,2, 
for some constant C and all sufficiently small > 0, and hence we get 

|X+(a+2(t2)) -X+(a+2(ii))| < |^(t2)+^(ti)| +C(V^+ V^)- (29) 

On the other hand, thanks to cr+2(ii) > cr* + e > cr* and the properties of X+, cf. Remark 2, we 
have 

|a+2(t2) -C7+2(tl)| < Co|X+((7+(t2)) -X+((7+2(tl))|, 

and combining this (29) gives 

|a+2(t2)-a+2(ii)| <Co|£(t2) + ^(ti)|+c(v/CM+v/C(i2))- (30) 

In the general case we introduce again two times ti and t2 by using (28) with (1+2 instead of ao, and 
argue as above. Moreover, the estimate 

|a_2(i2)-a_2(ti)| <Co|£(t2) + ^(ii)|+c(yCM+yCM)- (31) 
can be proven similarly. 
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Step 3: By construction, we have 

|(7_i(t2) - (7-l(tl)| < 2£, \(T+l{t2) - (T+l{h)\ < 2e, 

and Assumption 5 implies 



\iit2) - i{ti)\ < [ sup m ]\t2-ti\<Co\t2-ti\. 

\te[ti,t2] J 

The desired result now follows from the algebraic relation (26) as well as the estimates (27), (30), 
and (31). □ 

4.3 Compactness results and convergence to limit model 

In this section we finally pass to the limit v ^ and verify the validity of the limit model. We 
therefore write 

instead of r, Qi^ instead of g, Ui, instead of a, fnj^v instead of rrij, C,e,v instead of C,, 
and define the phase fraction by /Uj/ := m-i^^n — m-^i,. 

Theorem 29 (convergence to limit model along subsequences). There exists a sequence {i'n)n&n 
— >■ as n — >■ 00 as well as two Lipschitz functions ctq, /xq G C°'''^([0, T]) such that 



I II n— >oo „ 
- Cro\\c(I) > U, 



I II „ 

iMl'n - /^0||C(/) > U, 



(32) 



for each compact interval I C (0, T]. Moreover, we have 

n— >-oo 1 - fJ,o{t) 



QuAt, x) 



weakly* for all t > 0, and the triple {£, ctq, /^o) is a solution to the limit model in the sense of 
Definition 10. 

Proof. Convergence of a: We choose a sequence (SrOrteN with < < e* for all n and e„ — as 
n — )• 00. According to Lemma 27 and Lemma 28, there exist - for any given n - positive constant 
Co, Cn, Cin, and /3„ < 1 such that 



^u{to) < r 



sup CenA*) ^ C'n^^ 

te[to,T] 



and 



|(T^(i2)-a4tl)| <Co(|t2- til +C„(t^"+ sup J CenA^)) 

^ ^ ^ tG[tl,t2] ^ ^ 

holds for all n, all times to,ti,t2 G (0, T], and all sufficiently small > 0, where Cq is in fact 
independent of n. Moreover, making Co larger (if necessary) we can also assume that 



Cor^> / V^{t) 
Ju 



dt > T,f " 



{te[t,,T] : V,{t)>4-] 



holds for all > and n G N, and hence there exists for any choice of v and n a time 

Sn,v G [i*, + Co4^-''")] with V^{Sn,u) < T^"- 

For each n we next choose f„ > such that 
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In particular, using the abbreviations a^, := a^^, rrij^n ■= ruj^i,^, Cn := Cen.^-n) and Sn := Sn,u„ we 
have 

sup Cn{t) < < £n (33) 

t&[S„,T] 

as well as 

\(Tn{t2)-(Tn{h)\<Co\t2-h\ + {Co + l)en for ah h,t2 & [Sn, T] . (34) 

Let ^0 > be fixed and notice that Sn < to for almost all n. The compactness result from Appendix 
C, Proposition 35, guarantees the existence of a continuous function ctq defined on [to, T] and a 
not relabeled subsequence such that ||(7n — cro||c([to,T]) ^ as n — )■ oo. Moreover, by the usual 
diagonal argument we can extract a further subsequence such that — o"o||c(/) ~^ foi' any compact 
I C (0, T], and the estimate (34) implies that ctq is Lipschitz continuous on the whole interval [0, T]. 

Convergence of fj, and g: In what follows, we denote by Co any generic constant independent of 
n, and assume (without saying so explicitly) that n is sufficiently large. We also define 

a; := |(o-* + cj#), a := |((7* + a*) , 

and introduce a function fiQ as follows: For ao{t) G (— oo, a} we set Ho{t) = —1, and since we have 
Eje{-,o,+} ^j,nit) = 1 as weh as 

mo,n{t) + m+,n{t) < Cn{t) for en<a-a* 

we find 

\lJin{t)- ^Xoit)\<Cn{t). (35) 

Similarly, for ao{t) G [a, +oo) we set iJLo{t) = +1 and find again (35). In the case ao{t) G (a, a), we 
employ Lemma 20 - applied with e = imn{a* — a,a — a*}, which does not depend on n - to find 

m - ^-^^f^X.{anit)) - l±M)x+(a„(i))| < CoUt), 

and hence £{t) G [X_(cr*), X+(cr*)]. In particular, we can define Ho{t) G [—1, +1] as the unique 
solution to 

m = '-^^X.{ao{t)) - '-^f^X^iaom (36) 

and using the properties of X±, see Remark 2, we prove that (35) holds also in this case. 

In summary, we have now defined Hoit) for all t G [0, T], and (35) combined with (33) and 
Sn ^ yields ||/in — /^o||c(/) ^ as n — >■ 0. Moreover, the claimed weak* convergence of Qn is a 
direct consequence of ^n{t) + '^o,n(i) < £n ^ 0, see Remark 21. 

Verification of limit dynamics: Using Lemma 20 once more we find 

X_ (c7(t)) = £(t) for ao{t)<a#, X+{a{t)) = i{t) for ao{t)>a*. 
and this implies 

^io{t) = -1 for ao{t) < (7#, fio{t) = +1 for ao{t) > a*. (37) 
Combing these results with (36) we readily verify the algebraic relations 

{e{t), aoit), i^o{t)) G n, C{£{t), ao{t), /xo(t)) = 0, (38) 
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where Q and C are defined in (10) + (11). Therefore, and thanks to the properties of i and the 
functions X±, see Definition 5 and Remark 2, the pointwise identities (38) imply that fiQ belongs in 
fact to C°'^([0, T]) and satisfies 



a 



for almost all t € [0, T]. 

It remains to establish the dynamical relations from Definition 10. We first observe that Lemma 
25 yields /t = for almost all t with a{t) G (c7#, cr*) , and combining this with (37) we conclude that 
a{t) ^ {<7#, <T*} implies fi{t) = 0. Now let t be a time such that a{t) = o"^ and yu(t) G (—1, +!)• 
The set constraint (n, o") G H then implies &{t) > 0, and in the case of &{t) = we can employ 
Lemma 26 to show that £{t) < and hence /i(t) < 0. The derivation of /i(i) > for a{t) = a* is 
similar. □ 

Notice that Theorem 29 neither implies cr^^{0) — >• (To(0) nor /ii/„(0) — )• /Lto(O). This is not surprising 
because we expect, as explained within §2, that each solution with generic initial data exhibits a 
small initial transition layer. More precisely, if the mass at time t = is not yet concentrated in 
two narrow peaks, the systems undergoes a fast initial relaxation process during which a and fi may 
change rapidly. After this process, that means at some time of order at most tI~^ , < /? < 1, the 
dissipation is of order and our peak stability estimates imply that afterwards the state g,y can be 
described by two narrow peaks, which in turn are either transported by the dynamical constraint or 
exchange mass by a Kramers-type phase transition. 

The above arguments reveal that the limit functions ao and can (and in general they do) depend 
on the subsequence, or equivalently, on the microscopic details of the initial data. For well-prepared 
initial data, however, we can improve our result as follows. 

Theorem 30 (convergence for well-prepared initial data). For well-prepared initial data in the sense 
of Definition 8, we can choose I = [0, T] in (32). In particular, the whole family ((£, ctj/, Mi/))i^>o 
converges as u ^ Q to a solution to the limit model. 

Proof. By assumption, there exist values cTini G M and /iini G [—1, 1] such that a^{jd) — cjini as well 
as ^i/(0) — /Uini as — 0. Now let (((J„, ^n))neN sequence as provided by Theorem 29. Since 

the initial data are well-prepared, we can choose = in the proof of Theorem 29, see also Remark 
9. This implies 

II II , II II '^-^■O n 
iFn - Cro||c([0,T]) + ll/^n - /^0||C([0,T]) > U 

and hence cr„(0) — > (Jini as well as /x„(0) /Xmi. Since the limit model has only one solution with 

initial data (dmi, /^ini), see Proposition 34, we conclude that each sequence from Theorem 29 has the 
same limit, and standard arguments (compactness+uniqueness of accumulation points=convergence) 
provide the claimed convergence. □ 



A Solutions to the nonlocal Fokker-Planck equation 

In this appendix we show that the initial value problem to the nonlocal Fokker-Planck equation 
(FPi) and (FP2) is well-posed with state space 

P^(M) := jprobability measures on M with bounded variancej. 

We emphasize that all results derived in this section apply to arbitrary (i.e., uncoupled) parameters 
u>Q and r > 0. 

Our existence and uniqueness proof is based on a fixed point argument that allows to construct 
solutions to the nonlocal problem by iterating the solution operator of a linear PDE with a nonlinear 
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integral operator. The key ideas are as fohows. Let r, v be fixed and ^^ini be given. For any 
(7 G C( [0, T] ), we denote by 7^[o"] the solution to the linear PDE (FPi). In other words, for each a 
the function 7?.[cr] satisfies the initial value problem 

r5t7^[a](^, X) = i^^dln[a]{t, x) + d^(^{H'{x) - a{t))n[a]{t, x)), n[a]{0, x) = gini{x) (39) 

with a; G M and t G [0, T]. Using TZ, we now observe that the dynamical constraint (FPg) is equivalent 
to the fixed point equation a = S[a], where the operator S is defined by 

S[a]{t) := / H'{x)n[a]{t, x) dx + Ti{t). 

Notice that (FP2) implies (FP2) if and only if the initial data are admissible in the sense of 
J^XQini{x)dx = i{0). 

Our first result in this section employs Banach's Fixed Point Theorem in order to show that S admits 
a unique fixed point in the space of continuous functions. Afterwards we derive some bounds for 
these solutions which are uniform with respect to r and u. 

Proposition 31 (Existence and uniqueness of solution). For any t > 0, u > and all initial data 
Qini € P^(M) there exists a unique solution to (FPi)-/-(FP2). In particular, g is smooth in (0, T] xM as 
well as continuous in t with respect to the weakk topology in P^(M), and a is continuously differentiable 
on [0, T]. 

Proof. Operators and moment balances: For given o" G C( [0, 1] ), the existence, uniqueness and reg- 
ularity of 7^[u] can be established by adapting standard methods. For instance, [Fri75, Section 6, 
Corollary 4.2 and Theorem 4.5] guarantees the existence and uniqueness of smooth solutions under 
slightly stronger assumptions (boundedness of H'). For linearly increasing H', wc arc only aware 
of results concerning the stochastic Langevin equation rdx = (o"(t) — H'{x)^ dt + v^dW, sec e.g. 
[Fri75, Section 5, Theorem 1.1]. The solution TZ[a] to (39) is then provided by the resulting proba- 
bility distribution function for finding a particle at {t, x). We also refer to [JK098, ASZ09], which 
study the existence and uniqueness problem for similar equations in the framework of Wasserstein 
gradient flows. 

Using the PDE as well as integration by parts we verify the moment balance 

r-^ [ il;{x)Q{t, x)dx = [ ^p" {x)Q(t, x) dx + [ tp'{x){a{t) - H'{x))g{t, x)dx (40) 

for any tp with lV-'(-2^)| + |V'"(3^)I < ^(l -|- x^) and |'i/''(a;)| < C(l + \x\) for all x G M, and this implies 
the desired continuity of moments with respect to t. For tp(x) = 1 we obtain J Q{t, x)dx = 1 and 
with 'ip{x) = x^ we verify that 

/" x^Q{t, a;)dx < ( 1 -|- [ x"^ Qini{x) dx \ exp ( C^*^ + ll^lloo A ^ 

Jr \ Jr J \ t J 

Moreover, the choice ip = H'{x) reveals that the operator S is well defined since H'{x) grows linearly 
as a; — >■ ±00, see Assumption 1. 

Lipschitz estimates: We next consider two functions cri, (T2 G C( [0, T] ) , abbreviate Qi := TZ[ai], 
and introduce functions Ri and R2 by 

/X 
Qi{t, y)dy. 
-00 

The function R := R2 — Ri then satisfies 

TdtRit, x) = u'^dlR{t, x) + {H'{x) - (j2{t))d:cR{t, x) - {a2{t) - ai{t))gi{t, x). 
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In view of gi{t, •) G P^(M) we readily verify that 



x'^ \R{t, x)\ 

and this imphes R{t, •) G L^(M) for all t as well as 



a;->±oo 



> 0, 



S[a2]{t) - S[ai]{t) 



H'(x)d^R{t, x)dx 



[ H"{x)R{t, x)dx <C j \Rit, x)\ dx. 



In order to derive L^-bounds for R, we fix some e > and approximate the modulus function by 
he{r) := Vs + r'^. Thanks to -1 < h'^{r) < 1 and h'^{r) > for all r G M, we obtain the moment 
estimate 

T-^ / he{R{t, x)) dx<- [ H"{x)h,{R{t, x)) dx - ((72(t) - (7i(t)) / K{Rit, x))gi{t, x) 
ci* Jr Jr Jr 

<C I he{R{t, x)) dx + a2{t) - ai{t) , 
Jr 

where C := ||i?"||oo. Using the comparison principle for ODEs and passing to the limit e — >■ we 
therefore get 



/ \R{t, X)\dx <T~'^e^p{CT~H) [ C72(s)-C7l( 

Jr Jo 



ds, 



where we used that R{0, •) = holds by construction. 

Fixed point argument: The estimates derived so far ensure that 



S[a2]{t) - S[ai]{t) = C f a2{s) - (7i(s) 

Jo 



ds, 



and this implies that 5 is a contraction with respect to llfjUr '■= ^'^Vte[o,T] ^^P (— 2Ct) |cr(t)|, which 
is equivalent to the standard norm. The existence of a unique fixed point is therefore granted 
by Banach's Contraction Principle. Now suppose that S[o\ = a. Prom (40) with ip{x) = H'{x) 
and ip{x) = X we then conclude that a is continuously differentiable and that (PP2) is satisfied, 
respectively. □ 

Proposition 32 (Uniform bounds for solutions). Let t and v he fixed with < r < r and < i' < i'. 
Then, each solution to the nonlocal Fokker-Planck equation (PPi)-/-(PP2) satisfies 



sup ( \a 
te[o,T] 



{t)\+ [ X^Q{t, 

Jr 



as well as 



sup \\Q{t, 



x)dx\ <C 



C 



00 < ^ 



and 



f 



V{t) dt < Ct 



where C is some constant which is independent of r and u hut depends on H , t, v, I, and 

fj^X Qiyii(^X^dx. 
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Proof. Moment estimates: Due to the constraint (FP2), the moment balance (40) with ^(x) = x 
imphes 

T-^ [ x^Q{t, a;)dx < 2z/^ + 2||(7||oo|K||oo + 2c- c /" x^Q{t, x) dx, 

where c is chosen such that xH'{x) > c{x^ — l) holds for all x. Employing the comparison principle 
for scalar ODEs we therefore find 

/ x^Q{t, x)da; < 2i/2 + 2||(7||oo||^||oo + 2c+ / x^Qira{x)dx 
Jr Jr 

<C(l + ||a||oo). 
Moreover, by applying Holder's inequality to (FP2) we get 

1/2 / /• \ 1/2 



\a{t)\ < T\i{t)\ + (^J \H'{x)f Q{t,x)dx^ (^J g{t,x)dx^ 

/ x'^g{t, x) I 
\Jr 



1/2 

<C + C{ I x''p(t, x)dx 



where C is some constant independent of r and z^. The combination of both estimates gives 

||o'||oo < C^/l+ ||cr||oo, 

and the desired moment bounds follow immediately. 

L°° -estimate after waiting time Parabolic regularity theory implies that \\Q{t, •)||oo is well- 
defined for all t > 0, and thus we only have to understand how this quantity depends on t, r, v, and 
the initial data. To this end we fix to with <to <T, consider the function 

Mto(t) := sup \\VsQ{to + s, ■)\\oo, 

0<s<t 

and denote by C any generic constant that is independent of r, u and to- Using the rescaled heat 
kernel 



Kit, X) :=y^exp(- 



as well as Duhamel's Principle, any solution to (FPi)-|-(FP2) can be written as 

Q{to + t,x) = h,to{t, x) + h,to{t, x), 

where 

A,to(i, x) := / K{t, X - y)g{to, y) dy 
Jr 

and 

l2,to{t,x):=- [ [ K^{t-s,x-y)f{to + s,y)dyds, f{t,x):={H'{x)-a{t))g{t,x). 

Jo Jr 

The first term can be estimated by 

\h,to{t, x)\ < \\K{t, OIloo [ Q{to, y)dy < -.^ 

JR f \ I 

whereas for the second term we employ Holder's inequality to find 

\l2,toit,x)\<^J^(^J^K,it-s,y)^dy^ ' (^jjito + s,yfdy^ ' . 
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By direct computations we verify 

and using \H'{x)\ < C(l + Q{t, x)dx = 1 as well as the uniform moment bounds we get 

/ /(to + s, y)'dy<C||^(to + s, •)lloo(|(T(to + s)|' + l+ / y'^(to + s, y)dy\ 

JR \ JR / 

<Cs-^'^Mt,{s). 
The latter three estimates imply 

where we used Jq {t — s)~^^^ s"^/'^ ds = (1 — s)~^^^s~-^/^ ds < oo and that Mf^ is an increasing 
function in t. We therefore get 

/.II Ml ^ cv^ c^ym^ 

Vt\\g{to + t, •)l|oo<— + ^3/,^i/4 , 
and since an analogous estimate holds for all < s < t, we arrive at the estimate 

Mt„(i)<— + — ^7^-^. 

This implies 

Vt|k(to + i, •)lloo<Mto(t)<Cmax|^, (41) 

and for t = u'^t we get 

||£>(to + z^^r, •)||oo < ^• 

The claimed L°° -estimate now follows since to was arbitrary and C independent of to- 
Bound for dissipation: The energy balance (3) implies 

J V{t)dt = T(^£{iy^T) -£{T) + ait)i{t)d?j 
<t{S{u\) -£{T) + C), 
and from the definition of the energy (1), the above L°°-bounds, and H{x) < C{l + x^) we infer that 

£{v^t)<v^ I Q[v^T,x)hiQ[v^T,x)dx + C I [1 + x^)q[v^t, x)dx 

JR JR 

< In -5 I + C < C. 



In order to derive a lower for S{T), we assume (without loss of generality) that the global minimum 
of H is normalized to 0. The properties of H, see Assumption 1, then guarantee the existence of 
constants c > as well as x_ < and x+ > such that 

^ ^ - \ {x- x+y for x>0, 
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and hence we estimate 

J 7o(x) dx < J exp — j dx + J exp ^-'^ — j '^^ - 

where 'yo{x) := exp (— This implies 

S{T) = 1^^ f git, x) In ( ^^h^) dx 
Jr V 70 (a;) y 

V V 70 (a;) J Q{t, x) J Jr 

where we used that \az + l/ z> 1 holds for all z > 0, and the desired L^-estimate for the dissipation 
follows immediately. □ 

Remark 33. For initial data ^ni G L°^(M) we have 

sup II ^(t, OIloo < I V{t)dt<CT 

te[o,T] 2^ Jo 

for some constant C which depends only on H, f, v, I, J-^x"^ q\^\{x) dx , and ^'^||£»im||oo- 
Proof. In this case we can estimate 

A,o(i, a;) < llftnilloo / K{t, x)dx = ||ani||oo- 

Moreover, for < s < f < v^t we infer from (41) that 

V^\\q{s, OIIoo < Mo(s) < Mo{t) < ^ 

and this implies 



^2,o(i, x) < -^^^ 1^ (t - s)-'/y\\g{s, Ollood. < ^ 



C 

2 ■ 



'0 

The claimed L°°-estimate now follows from summing both inequalities (for < t < and using 
Proposition 32 (for i/^r < t < T). Moreover, the L^-bound for the dissipation can be derived as in 
the proof of Proposition 32. □ 

B Solutions to the limit model 

We prove that the initial value problem for the limit model has always a unique solution. 

Proposition 34. For any £ as in Assumption 5, and any given initial data a(0) and //(O) with 
(£(0), (t(0), ju(0)) G ft, there exist two functions a and v on [0, T] such that 

1. both a and fu, are continuous and piecewise continuously difjerentiahle, 

2. both functions attain the initial data, 

3. the triple {£, a, ji) is a solution to the limit model in the sense of Definition 10. 
Moreover, u and v are uniquely determined by £, cr{0), and f^{0). 
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Figure 12: Cartoon of the piecewise smooth vector field V+ (arrows) on the set S (gray area) as used in the 
proof of Proposition 34. For given initial data from S, there exists a unique integral curve which is continuous 
and piecewise continuous differentiable. 



Proof. We observe that 



{I, a, n) eO. 



(m, a) e 



where the closed set S is defined by 
S := {—1} X ^— oo, a 



# 



U (-1, +1) X 



# 



U {+1} X [a#, +oo) , 



see Figure 12 for an illustration. Moreover, for each point (^, a) e E there exists a unique value for 

£ such that C{i, a, //) = 0. Wc proceed with discussing three special cases: If i{t) = £{0) holds for 
all t G [0, T], then the unique solution to the limit model is given by a{t) = cr{0) and fi{t) = /u(0). 
In the case of i{t) > for all t G (0, T), we argue as follows. By reparametrization of time, we can 
assume that £{t) = 1. The pointwise constraint C(^£{t), a(t), //(i)) = then implies that any solution 
to the limit model satisfies 

for almost all t G [0, T], where the vector field V+ : S ^ is defined by 
(x+(a)-X_(a))~\ 



V+(/x, a) 



for 



0, 



1 < /X < +1 and a = 
for all other points in S. 



Since V+ is piecewise continuously differentiable with derivative on S, there exists a unique contin- 
uous integral curve emanating from the initial data, and this integral curve is obviously piecewise 
continuously differentiable. The arguments for the third case, that is £{t) < for all t G (0, T), are 
entirely similar. For arbitrary £, we introduce times = Tq < Ti < ... < T/v = T such that for any 
i = 1...N and all t G (Tj-i, Ti) we have either £(t) < 0, or i{t) = 0, or £{t) > 0. The assertion now 
follows by iterating the arguments for the special cases. □ 



C Non-standard compactness criterion for continuous functions 

In the proof of Theorem 29 we utilize the following, non-standard compactness result in the space of 
continuous functions. 

Proposition 35. Let I be some compact interval, g G C(/) a continuous function on I, and (cn)„gi!j 
be a positive sequence with c„ — >■ as n — >■ 0. Moreover, let {fn)neN ^ ^(-^) ^ bounded sequence 
such that 

\fn{t2) - fnitl)\ < \9{t2) - gitl)\ + C„ (42) 
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holds for all ti,t2 G / and all n G N. Then, the sequence {fn)nen *^ compact in C(/) and hence 
equicontinuous. 

Proof. Without loss of generality we assume that / = [0, 1]. Since the sequence (/n(t))„gf^ C M 
is compact for any i € M, we can - by the usual diagonal argument - extract a (not relabeled) 
subsequence, such that /n(i) converges as n ^ oo for alH G / fl Q. Our assumptions imply that the 
function : iCiQ^M. defined by 

/oo(i) := lim for all telnQ, 

n— >-oo 

satisfies 

|/oo(i2) - /oo(ti)| < 15(^2) - 9iti)\ for all ti,t2 € I n Q, 

and we conclude that foo admits a unique continuous extension f^o G C(/), which obviously satisfies 

\foo{t2)-foo{ti)\<\g{t2)-9{ti)\ for all ti,t2el. (43) 

We next show that fn converges to /oo as ra ^ oo strongly in C(/). To this end let 5 > be fixed. 
Exploiting the continuity of g as well as (42) and (43) , we first choose no G N and iV G N such that 

n>no, \t2-ti\<j^ =^ \fn{t2)- fn{tl)\ + \foo{t2)- fooitl)\<S/2, (44) 

where both no and N can depend on S. We next divide I = [0, 1] into N subintervals of length 
that means we introduce 

= to < ti < t2 < ... < tN = i with tj := j/N. 

For each t e I there exists j = j{6, t) G {0, 1, ...N} such that \t — tj\ < 1/iV, and (44) ensures that 

n > no =^ \fn{t) - fn{tj)\ + \foo{t) - foo{tj)\ < 6/2. 

We finally choose ni such that 

n>ni ^ sup \fn{tj) - foo{tj)\ < S/2, 

ie{o,i,...iV} 

and combining the latter two implications gives 

n>max{no,ni} =^ \fn{t) - foo{t)\ < 6 

for all t & I. Since 5 was arbitrary, we have thus proven that ||/n — /oo||oo ^ as n — >^ oo, and the 
equicontinuity follows from the Arzela-Ascoli Theorem (e.g. [DiB02, Proposition 19.1]). □ 
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